NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA 
SUBMARINE  HYDRAULIC  CONTROL  ANALYSIS. (U) 
JUN  BO  M  J  BOWER 


ADA090091 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


DTIC 

ELECTE 

OCT  9  1980 


'/  / 


U- 


\ 


THESIS 


1 _ j 


o 

o 


/  SUBMARINE  HYDRAULIC 
r" 

CONTROL  ^ANALYSIS, 

by 

^  /  Michael  J .. 

Bower 

d/  • 

/-  Jun-®* 

sd 

Thesis  Advisor: 

T.  M.  Houlihan 

Approved  for  public  release;  distribution  unlimited 


D 

* 


■a 

i 


SECURITY  CLASSIFICATION  OF  THIS  RAPE  fAaw  Pin  hioN) 


1  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS  j 

BEFORE  COMPLETING  FORM  1 

4.  TITLC  (mi  SuMHta) 

SUBMARINE  HYDRAULIC  CONTROL  ANALYSIS 

Master  s  Thesis; 

June  1980 

RERFORMINO  ORB.  RERONT  NUMBER 

T.  AUTHOR**) 

Michael  J.  Bower 

■■■■ 

S.  RERFORMINO  ORGANIZATION  HAMS  ANO  ADDRESS 

Naval  Postgraduate  School 

Monterey,  CA  93940 

10.  RROORAM  ELEMENT,  RROlECT.  TAM 

AREA  A  MORN  UNIT  NUMBERS 

II.  CONTNOLLINO  OFFICE  MAMS  ANO  AOORESt 

Naval  Postgraduate  School 

Monterey,  CA  93940 

IS.  RBRORT  DATE 

June  1980 

65  _ 

It  uOM|ToRin6  AOEHCY  NANS  S  AOORESt<lf  HIImM  frmm  Cmlrmltht*  OWem) 

IS.  SECURITY  CLASS.  (•«  Ala  r*M«) 

Unclassified 

ISa.  dcCL ASSlFlCATION/OOMNaRADlNO 
SCHEDULE 

f«.  WfTHltuTlON  IT  ATKMtNT  ft  •  KfH) 

Approved  for  public  release;  distribution  unlimited 

IT.  DISTRIBUTION  STATEMENT  (••  »•  ataMal  «im4  In  BtM S  JS.  II  Hffmml  tnm  *N»IJ 

is.  surrlbmsntary  notis 

If.  K  tv  MONOS  fCmiliM  an  •14a  II  mMMr  «RM  limHfr  *r  Mm*  mMmI 

Hydraulics,  submarine  hydraulics,  controls,  controls  analysis 

SO.  ABSTRACT  fCmtl w»»  m  kwh  AM  II  mimmn  M  lAaaMfr  Sr  Mm*  mini 

/  -A  mathematical  model  was  developed  to  include  the  effects 
in  the  submarine  hydraulic  system  dynamic  performance  analysis. 

The  project  was  undertaken  in  an  effort  to  demonstrate  the 
necessity  of  coupling  the  entire  hydraulic  power  network  for 
an  accurate  analysis  of  any  of  the  subsystems  rather  than  the 
current  pr~  tice  of  treating  a  component  loop  as  an  isolated 
system.  It  was  intended  that  the  line  model*  could  be  coupled 

DO  ,'2Tr»  WTCy:  eoition of  i nov ••  ••  oMOLCTs  Unclassified 

S/N  0 103*014*  SA0I  I  SECURITY  CLASSIFICATION  OF  TNI*  FAOC  fVMN  brnm" SiiarM) 


Unclassified 


PXEC mil* ib  lai'iiut  i 


with  individual  elements  in  modular  computer  programs  to 
facilitate  design  and  alteration  decisions. 

The  line  model  was  developed  and  coupled  with  a 
submarine  driving- system  servo- loop  in  order  to  demonstrate 
the  above  premises.  However,  due  to  programming  difficulties 
additional  work  is  necessary  to  successfully  complete  the 
analysis . 

/ 


D0,SS3a  UT3 

S/N  0102-014-6601 


Unclassified 

IICu*i*v  Cb*Min«*ri«M  rm»  rntgBtmmi  hn 


2 


Approved  for  public  release;  distribution  unlisted 


Submarine  Hydraulic  Control  Analys 


as 


by 


T  j  Michael  J.  Bower 

sz 

S|“rewncf  fo?lthe1S^”n|f0f  the 
master  of  science  in  mechanics  engineering 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1980 


Author : 
Approved  by: 


Thesis  Advisor 


Second  Reader 


* *  CAM. 


ABSTRACT 


A  mathematical  model  was  developed  to  include  line 
effects  in  the  submarine  hydraulic  system  dynamic 
performance  analysis.  The  project  was  undertaken  in  an 
effort  to  demonstrate  the  necessity  of  coupling  the  entire 
hydraulic  power  network  for  an  accurate  analysis  of  any  of 
the  subsystems  rather  than  the  current  practice  of  treating 
a  component  loop  as  an  isolated  system.  It  was  intended 
that  the  line  model  could  be  coupled  with  individual  elements 
in  modular  computer  programs  to  facilitate  design  and  alter¬ 
ation  decisions. 

The  line  model  was  developed  and  coupled  with  a  submarine 
driving- system  servo- loop  in  order  to  demonstrate  the  above 
premises.  However,  due  to  programming  difficulties 
additional  work  is  necessary  to  successfully  complete  the 
analysis . 
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I.  INTRODUCTION 


A.  GENERAL  DISCUSSION 

The  vast  majority  of  hydraulic  systems  currently  in  use 
are  of  a  dedicated  type.  That  is,  they  are  designed  for  a 
specific  function  with  consistent  and  well  defined  constraints. 
In  most  cases  they  are  further  limited  to  a  loop  in  which 
the  pump  supplies  power  to  a  single  power  element.  Such 
systems  can  be  delivered  as  compact  units  with  a  minimum  of 
connecting  piping.  Some  examples  of  dedicated  systems  in 
Naval  use  are:  surface  ship  steering  units,  gun  mount 
train  and  elevation  units,  and  ammunition  conveyor  systems. 
Although  submarines  have  employed  dedicated  systems  in 
the  past,  all  modem  U.S.  submarines  use  constant  pressure 
systems  [Ref.  1,  Tab.  561-1].  The  constant  pressure  systems 
contrast  with  the  dedicated  type  in  that  the  supply  and 
return  headers  extend  the  entire  length  of  the  ship  and  serve 
a  multitude  of  power  elements.  The  total  piping  length  of 
these  systems  may  exceed  the  length  of  the  ship  by  a  factor 
of  four  or  more. 

B.  THE  NEED  FOR  A  LINE  MODEL 

For  the  packaged  configuration  of  the  dedicated  system, 
with  its  relatively  short  piping  lengths,  the  current  practice 
of  neglecting  line  dynamics  and  lumping  the  system  volume 
has  produced  acceptable  analytical  results.  However,  even 
in  systems  with  low  line-volume  to  total-volume  ratios,  the 
computed  values  of  the  hydraulic  natural  frequency  (w^) 


are  frequently  in  excess  of  40%  higher  than  measured  values 

[Ref.  2].  Since  co^  sets  the  upper  limit  on  system  response 
and  is  affected  by  the  volume  as  well  as  line  dynamics,  a 
model  which  includes  these  factors  would  undoubtedly  produce 
more  accurate  results . 

It  was  the  intent  of  this  project  to  demonstrate  that 
the  submarine  hydraulic  system,  with  its  extremely  long  lines, 
requires  that  the  line  dynamics  be  included  to  achieve  any 
degree  of  accuracy  in  analytical  solutions. 

C.  POTENTIAL  FOR  SYSTEM  MODEL 

Discussions  between  the  author  and  Mare  Island  Naval 
Shipyard  Hydraulic  Design  Division  personnel  (code  260)  > 
as  well  as  a  Naval  Sea  Systems  Command  representative 
[Refs.  2  and  3],  disclosed  that  there  is  currently  no  computer 
program  to  model  the  dynamic  behavior  of  a  submarine  hydraulic 
system.  Such  a  program  would  be  useful  in  the  design  of  new 
systems,  both  for  sizing  of  lines  and  elements  and  for 
determining  the  response  of  components  served.  For  existing 
systems,  it  would  allow  a  determination  of  the  effects  of 
alterations  prior  to  tneir  implementation. 

Additionally,  a  computer  program  of  a  dynamic  system 
model  would  be  advantageous  in  analyzing  vibrations  trans¬ 
mitted  through  the  system.  The  "noise"  generated  by 
vibrations  is  a  critical  element  of  the  "Submarine  Sound 
Quieting  Program,"  while  the  fatigue  due  to  vibrations  is 
important  to  the  life  of  mechanical  elements  in  the  system, 
especially  the  seals  [Refs.  5,6,  and  7], 
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The  development  of  a  dynamic  line  model  which  could  be 
programmed  in  a  modular  format  and  coupled  with  models  for 
the  individual  elements  would  facilitate  the  creation  of 
such  an  overall  system  model. 
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II .  SYSTEM  MODELED 


A.  MAGNITUDE 

It  was  recognized  early  in  the  undertaking  that  it  would 
not  be  possible  to  develop  a  model  for  an  entire  submarine 
hydraulic  system  in  the  time  available.  For  this  reason, 
the  stem  diving  planes  servo-loop  was  chosen  for  the 
primary  development  (Fig.  1).  This  sub-system  was  considered 
to  be  representative  of  one  of  the  more  complex  loops  in  the 
overall  system.  (  a  loop  being  defined  as  the  fluid  path 
through  a  control  valve  -  to  the  power  element  -  back 
through  the  control  valve) .  Successful  modeling  of  the 
stem  planes  servo-system  would  readily  permit  expansion 
of  the  system  to  include  the  fairwater  planes  and  rudder 
loops  in  parallel  operations.  With  these  three  sub-systems 
modeled,  the  remainder  of  the  hydraulic  system  could  be 
relatively  easily  developed. 

B .  REALISM 

An  operational  ship's  hydraulic  system  was  deemed 
essential  as  the  basis  for  the  model,  as  it  would  allow  the 
computed  parameters  to  be  compared  with  actual  performance 
data.  The  hydraulic  plant  of  USS  PERMIT  (SSN  594)  was 
chosen  as  the  system  to  be  modeled.  It  was  considered 
representative  of  a  modem  submarine  hydraulic  system  and 
the  technical  information  was  readily  available  at  Mare 
Island  Naval  Shipyard. 


ANALYZERUtraa-^  X  t  III  SERVOVALVE 
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STERN  PLANES  SYSTEM 


Ill .  HYDRAULIC  FLUID 


A.  TYPE 

The  hydraulic  oil  used  in  all  U.S.  submarine  internal 
high  pressure  (3000  psi)  systems  is  Navy  symbol  2190TEP. 

This  is  a  petroleum  based  fluid  conforming  to  military 
specification  MIL-L-17331  [Ref.  8].  It  exhibits 
approximately  the  same  characteristics  as  an  SAE-30-Detergent 
motor  oil. 


B .  PROPERTIES 
1.  Density 

The  density  (p)  of  the  oil  was  obtained  from  Ref.  9 
in  terms  of  specific  gravity  (SG) .  The  values  supplied 
were 

SG  -  0.896  at  74°F. 

SG  -  0.0901  at  60°F. 

The  density  of  water  at  60  degrees  Fahrenheit  is 
-5  3 

9.3475  x  10  (slug/in  )  and  is  assumed  to  be  the  reference 
for  specific  gravity  unless  otherwise  specified  [Ref.  10] . 
Specific  gravity  is  defined  as  the  ration  of  the  density  of 
the  oil  at  the  given  temperature  to  that  of  water  at  60  degrees 
Fahrenheit. 


SG 


p 


(1) 


Solving  for  the  density  of  oil  at  temperature  (T) 
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yields 


p  -  8.4221  xlO"5  slug/ in3  at  60°F 

p  -  8.3754  xlO'5  slug/ in3  at  74°F 

2.  Cubical  Expansion  Coefficient 

The  cubical  expansion  coefficient  (a)  is  the  fractional 
change  in  volume  of  the  fluid  due  to  a  change  in  temperature  and 
is  used  in  the  determination  of  density,  viscosity  (u,v), 
and  bulk  modulus  (6)  values. 

Density  values  arrived  at  in  the  previous  section 
were  used  to  obtain  a  value  for  the  cubical  expansion  coefficient. 

Using  the  form  given  in  Ref.  11, 

p(P,T)  -  pQ[l  +  (P-P0)/B  -  a(T-To)]  (2) 

where  the  subscript  (Q)  denotes  reference  values  of  one 
atmosphere  of  pressure  and  60°F  temperature. 

Since  the  values  for  the  density  of  the  oil  were 
obtained  at  one  atmosphere 

P  -  Po 

equation  (2)  reduces  to 

p  -  pQ  [1  +  a  (T  -  T0)]  (2a) 

Inserting  appropriate  values  into  (2a)  and  solving  for  a 
yields 

a  -  3.9638  x  10‘4  °F'1 

3.  Viscosity 

a.  Kinematic  Viscosity 

Reference  12  specifies  that  2190TEP  oil  shall 
have  a  kinematic  viscosity  (v)  of  190SSU  (Saybolt  Seconds 
Universal)  at  a  temperature  of  130°F. 


(3) 


Using  the  standard  conversion  [Ref.  14j  for 

kinematic  viscosity,  where  the  time  (t)  is  greater  than 

100  seconds  (SSU) 

v  -  0.220  t  -  (135/t) 

yields  _ 

v  -  41.09  cs  at  130°F 

v  -  92.03  cs  at  100°F 

Converting  from  centistokes  (cs)  to  English  units  was 

accomplished  with  the  conversion  factor  (Ref.  15)  of  1.549  xlO 
2 

(in  /sec)/cs.  Thus  the  values  obtained  for  kinematic  viscosity 
were 

v  -  0.06367  in2/sec  at  130°F 
v  -  0.1427  in2/sec  at  100°F 

b.  Absolute  Viscosity 

Knowing  values  of  density  and  kinematic  viscosity 

permits  the  determination  of  absolute  viscosity  (u) . 

U  -  pv  (4) 

Using  the  form  presented  in  Ref.  16, 

u  -  yQe  -  X  (T  -  Tq)  (5) 

at  one  atmosphere  of  pressure,  where  \  is  the  temperature 

coeffcient  of  absolute  viscosity  with  units  of  (°F~^) . 

Combining  equations  (4)  and  (5)  and  then 

solving  for  lambda  gives 

\  -  In  (P0v0/pv)/(T-T0)  (6) 

Arbitrarily  setting 

T  -  100°F  and  T  -  130°F 
o 

in  conjunction  with  values  of  density  and  kinematic 
viscosity  at  these  temperatures  yields 
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Tq  -  100  T  -  130°F 

pQ  -  8.2886  xlO”5  p  -  8.1884  xlO"5 

vQ  -  1.427  xlO-1  v  -  6.367  xlO"2 

Substituting  these  values  into  equation  (6)  produced  a 
value  of 

X  -  2.731  xl0“2  °F"1 

Reference  17  presents  the  pressure  conversion 
for  absolute  viscosity  as 

log10(y/yo)  -  cP  (7) 

However,  as  the  other  functions  are  related  by  the  difference 
in  pressure  divided  by  the  bulk  modulus  (AP/8) ,  it  would 
seem  reasonable  to  assume  that  c  in  equation  (7)  has  a  form 
of  (c^/ B) .  This  would  then  yield  an  equation  for  absolute 
viscosity  at  constant  temperature  of 

u  »  yQe  -(c-j^&P/B)  (8) 

Lacking  any  other  information,  c^  was  set  equal  to  one, 
and  a  final  form  for  absolute  viscosity  was  formulated  as 

V  -  uoe«P/8)  -  X(T-T0>)  (9) 

where 

P  -  P-po 

P  ■  reference  pressure  (o  psig) 

P  -  system  pressure  (psig) 

8  *  system  effective  bulk  modulus  (psi) 

T  -  system  temperature  (°F) 

T  ■  reference  temperature  (°F) 
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4 .  Bulk  Modulus 


Although  the  literature  on  hydraulic  systems  indicated 
that  the  hulk  modulus  (6)  would  be  supplied  with  the  oil, 
visual  inspection  of  shipping  containers  and  bills  of 
lading  proved  futile.  Likewise,  the  military  specifications 
(MIL-L-17331)  failed  to  produce  any  worthwhile  information 
in  this  regard.  A  nominal  value  of  bulk  modulus  for  the 
oil  was  therefore  assumed  to  be  220,000  psi  ,  as  recommended 
in  Ref.  18. 

C .  LIMITATIONS 

In  a  further  effort  to  maintain  the  scope  of  the  project 
within  reasonable  constraints,  it  was  decided  to  treat  the 
temperature  as  constant.  Inferences  drawn  from  the  Naval 
Ships'  Technical  Manuals  previously  cited  and  the  author's 
own  experience,  indicated  that  a  temperature  of  130°F  was 
a  reasonable  choice.  With  this  value  of  temperature,  the 
equations  for  the  other  properties  were  simplified. 

From  equation  (2) 

p  -  8.188  xlO"5  +  8.4221  xlO'5(P/S)  (10) 

and  equation  (9)  yields 

u  -  5.213  xlO"6  e(P/8)  •  (11) 
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IV.  MODEL  DEVELOPMENT 


A.  SYSTEM  OPERATION 

Submarine  steering  and  diving  control  surfaces  (rudder 
and  diving  planes)  are  positioned  through  mechanical  linkages 
by  hydraulic  rams.  Each  surface  has  normal  as  well  as 
alternate  sources  of  hydraulic  and  electrical  power  to  ensure 
reliability.  For  the  purpose  of  model  development,  only  the 
normal  mode  of  the  stem  planes  was  considered. 

The  normal  mode  of  operation  for  the  SSN  594  class 
of  submarine,  employs  a  position  control  electro-hydraulic 
servo-loop  (Fig.  1).  In  this  type  of  system  (Fig. 2),  control 
station  operator  manually  positions  a  stick  to  achieve  a 
desired  angle  of  inclination  of  the  planes.  This  produces 
an  electrical  signal  which  is  summed  with  the  feedback 
signal  from  the  position  output  to  create  an  error  signal. 

The  error  signal  is  amplified  and  transmitted  to  the  torque 
of  the  servovalve.  The  torque  motor  displaces  the  wand  spring 
of  the  servovalve-pilot,  which  causes  an  imbalance  in  pilot 
pressure,  which  in  turn  shifts  the  main  spool  allowing  the 
servovalve  to  pass  oil  to  position  the  ram.  As  the  desired 
position  is  attained,  the  error  is  reduced  and  the  pilot 
returns  the  spool  to  a  neutral  positions ,  blocking  the  flow 
of  oil. 
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Figure  2  Stem  Planes  Control  Circuit 
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B .  SERVO-SYSTEM 

The  area  of  interest  for  the  model  development  was  the 
hydraulic  fluid  flow  in  the  servo-loop  as  depicted  in  Figure  3. 
By  inputting  the  valve  spool  motion  and  ram  load  force  as 
explicitly  defined  functions,  the  influences  of  the  external 
factors  involved  in  the  feedback  loops  (the  servovalve  torque- 
motor  and  pilot  valve,  and  the  ram  loading)  were  minimized. 

This  permitted  the  generation  of  flow  and  pressure  equations 
around  the  loop  with  a  minimum  of  complications  due  to  the 
ancillary  characteristics. 

1 .  Servovalve 

The  two  servovalves  commonly  in  use  in  the  Submarine 
Force  are  both  manufactured  by  the  Sargent  Corporation  and 
are  designated  Sanders  Model  SV  438-10P  and  SV  438-15P 
[Ref.  19].  Design  and  performance  of  both  valves  are 
similar.  The  SV  438-10P  which  is  installed  in  the  (SSN-594) , 
was  used  for  this  simulation.  Figure  4  is  a  representative 
drawing  of  the  actual  valve  employed  -  a  four-land  four-port 
element.  However,  comparison  with  the  valve  depicted  in 
Figure  3  -  a  three-land  four-port  model  -  will  show  that 
the  flow  is  essentially  the  same  in  both  valves .  It  was 
considered  easier  to  follow  the  flow  paths  of  the  valve  in 
Figure  3. 

All  technical  data  for  the  servovalve  was  obtained 
from  Ref.  20. 

a.  Flow  Paths 

The  pilot  valve  flow  was  disregarded, as  it  has 
its  own  supply  and  return  lines . 
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the  return 


Designating  the  supply  flow  as  Q^n, 
flow  as  Qr,  and  the  other  flows  by  their  respective  nodal 
points  (Fig. 3),  the  flow  equations  were  derived.  Using 


the  form  given  in  Ref.  21, 

EQin  -  E(W  =  dVdt  +  <V8>dP/dt  (12) 

the  results  were : 

Qin  "  Ql  +  Q2  (13) 

Qr  “  ^3  +  Q4  (14) 

Q1  =  Q3  +  Q5  +  Vl/31  +  h  (15) 

Q2  =  Q4  +  Q6  +  v2^2/32  +  ^2  •  (16) 


Here  the  volume  across  the  lands  of  the  valve  was  approximated 
as  zero.  The  superscript  (*)  represents  differentiation 
with  respect  to  time. 

Since  the  valve  is  manufactured  from  high 
strength  alloys  , 


as  shown  in  Appendix  A. 

In  addition  to  summing  the  flows  through  the 
valve,  the  flow  rate  was  also  defined  through  the  lands  for 
three  different  valve  positions. 

(1)  Null  Positions.  In  the  centered  (null) 
position,  a  critical-center  servovalve  will  have  leakage 
past  the  lands.1  This  flow  is  described  in  terms  of  null 


^or  a  comprehensive  coverage  of  all  terms  used  in  this 
paper,  the  reader  is  referred  to  Ref.  2. 
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(17) 


flow  pressure  coefficient  (KCQ) ,  as 

Q  -  %Kco(AP). 

The  four  node  flows  for  Xv  *  0  were 


*  *Kca<Pin-p> 

(18) 

«2  *  *Kco<Pin-p2> 

(19) 

^3  '  «eo<pl-pr> 

(20) 

«4  *  *Kco<P2-pr> 

(21) 

with 

P3  "  P4  "  Pr- 

See  appendix  B  for  values  Kcq  used  for  the  model. 

(2)  Valve  Spool  Off-Centered  Left .  The  valve 
travel  (Xv)  was  arbitrarily  defined  as  positive  to  the  left 
(Fig. 3).  Employing  the  orifice  flow  equation  [Ref.  22  ]  for 
the  fluid  passage  through  nodes  one  and  four,  and  an  empirical 
relationship  for  nodes  two  and  three,  the  flow  equations  were 


determined  to  be 

Qx  -  cdA/  2<"pin'Pl)/pl  (22) 

Q2  -  V*<W2  <ptn-p2>  <23> 

Q3  *  %Kcoe-X(Xv/Xvm)2  (PrPr)  (24) 

Q4  -  CdA/2(p2-Pr)/Pr  (25) 

where  Cd  =  0.61  for  full  periphery  rectangular  ports  and 

A  -  ird1Xv  (26) 

with  d^  defined  as  the  land  diameter  of  the  spool  valve. 

Since  a  value  of  leakage  flow  was  not 


determinable  from  Ref.  20,  the  empirical  relationships  of 
Appendices  C  and  D  were  developed  using  Ref.  23  as  a  guide. 
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(3)  Valve  Spool  Off-Centered  Right.  The 
equations  for  the  "negative"  flow  direction  were  similar  to 
those  developed  in  the  previous  section,  with  the  obvious 
exception  that  Q2  and  became  the  principal  flows.  In 
order  to  maintain  their  flows  as  positive  values ,  Xv  was 
squared  and  included  in  the  radical. 

The  form  for  Xv  less  than  zero  was 


Ql  -  V‘W2  <pln-pl>  <27> 

Q2  "  cdM  /2x£(pln'p2)/p2  <28) 

Q3  -  Cdw/  2X2(P1-Pr)/o2  (29) 

«4  *  <3°> 

with 

W  =  ^  (31) 


and  all  other  variables  as  defined  in  previous  sections. 

2 .  Lines 

The  line  equations  were  generated  with  the  intention 
of  accounting  for  all  factors  influencing  the  pressure  and 
flow  in  the  system.  The  schematic  diagram  of  the  actual 
system  employed  in  this  development  is  shown  in  Figure  5. 

For  development  purposes ,  only  the  odd  numbered  nodes  are 
shown  here  as  the  even  numbered  side  equations  were  identical, 
a.  Flow  And  Pressure 

Disregarding  all  but  the  main  line  piping,  the 
flow  path  from  the  servovalve  to  the  ram  was  as  depicted  in 
Figures  3  and  5  were  identified  by  the  same  numbers  .  The 
dynamics  of  these  segments  were  then  defined  at  each  point 
when  viewed  as  in  Fig.  6. 
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FIGURE  6  LINE  NODE  IDENTIFICATION 

Point  A  corresponds  to  nodal  point  (5)  of  Figure  3 
and  is  considered  to  be  located  at  the  entrance  of  the  pipe 
at  its  junction  with  the  servovalve.  Similarly,  point  B  is 
located  at  the  connection  of  the  two  different  diameter  pipes, 
in  the  entrance  of  the  downstream  pipe.  It  is  equivalent 
to  node  (7)  of  Figure  3.  Location  C  is  just  beyond  the  end 
of  the  line  segment  at  the  entrance  to  the  cylinder, 
identified  as  node  (9)  in  Figure  3. 

Utilizing  these  definitions  for  the  nodal  points 
and  viewing  the  direction  of  flow  as  always  being  from  A  to 
C,  the  dynamic  equations  for  each  element  were  developed. 

(1)  Node  Five .  The  pressure  and  flow  equations 
for  node  five  were  derived  from  Ref.  24. 

Q5  -  Q?  +  V57f>/B57  +  V57  (32) 

p5  -  pi  -  P5K50^/2A^  (33) 

The  double-subscripted  variables  were  employed 
to  denote  average  values  through  the  line  segment,  while  the 
single-subscripted  variables  represent  values  at  the  particular 
node .  Thus 
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A5L57 


~  1 

(34)  i 


with  being  the  actual  length  of  pipe  from  node  (5) 
to  node  (7)  and 

A5  -  ttd|/4  .  (35) 

As  can  be  seen  in  Figure  (6) ,  the  diameter  of  segment  A-B 
(5-7)  is  constant  and  equal  to  the  diameter  at  (5) .  The 
total  volume  and  actual  length  will  be  developed  further 
in  later  sections  of  this  report. 

l 

The  average  pressure  in  the  interval  was  j 

f 

P57  *  (P5  +  P7)/2  (36) 

with 

P57  *  ^57dC  '  <37> 

The  change  in  volume  with  respect  to  time  (V)  I 

was  set  equal  to  zero  for  all  line  segments ,  because  the 
lines  are  made  from  high  strength  materials  which  do  not 
permit  any  significant  expansion  or  contraction.  The 
exception  to  this  would  be  in  cases  where  there  are  flexible 
sections  installed  (  Aeroquip  -type  hoses) .  There  are 
numerous  locations  in  the  system  where  these  are  installed, 
but  not  in  the  servo-loop  modeled.  In  the  instances  where 
flexible  sections  are  incorporated  in  the  system,  an  expression 
for  the  change  in  volume  would  have  to  be  developed  based  on 
the  properties  of  the  material  used. 

The  value  of  Kg  was  based  on  an  abrupt  change 
in  the  cross-sectional  area  of  a  pipe, and  the  development 
is  shown  in  Appendix  E.  As  used  in  the  above  equations,  the 
sign  of  Kg  changes  with  the  direction  of  flow  such  that  the 
second  term  of  equation  (33)  will  always  reflect  a  pressure 
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loss  for  the  referenced  point.  Thus,  when  flow  is  positive, 
as  shown  in  Figure  (6),  P1  will  be  greater  than  P^,  but  when 
the  flow  is  reversed  (negative  in  Figure  (6))  Pg  will  be 
greater  than  P^. 

Qy  is  the  flow  at  node  (7)  and  will  be  defined 
in  a  subsequent  section.  P^  is  the  pressure  in  the  main 
servovalve  cavity  and  is 

Px  -  /Vt  (38) 

where  was  introduced  in  section  IV.B.l.a. 

(2)  Node  Seven.  The  flow  through  node  seven  is 
exactly  the  same  as  that  through  node  five,  with  subscripts 
changed  to  reflect  the  appropriate  line  segment. 

Qy  =  Qq  ^79^*79/^79  ^79  (39) 

P79  *  /P 79dt- 

The  pressure  at  point  (7),  on  the  other  hand, 
must  reflect  the  losses  due  to  flow  in  line  segment  (57) . 

This  was  accomplished  with  a  combination  of  the  equations 
in  Ref.  24  with  those  in  Appendix  E: 

P -j  *  Fj  “  32vj2^Le^yQgy/AjD^  -  1 . 14p^yQ^y/Aj  -pyKyQy/2A^ 

The  second  term  on  the  right  side  of  equation 
(40)  is  the  Hagen-Poiseuille  term  for  fully  developed  laminar 
flow  in  pipes.  The  third  term  accounts  for  developing  laminar 
flow.  The  final  term  accounts  for  the  entry/exit  loss  at 
cross-sectional  area  changes,  as  shown  in  Appendix  E. 

All  three  terms  must  reflect  a  loss  of 
pressure  for  flow  in  either  direction.  Since  Q57  assumes 
a  negative  value  for  flow  reversal,  it  presents  no  problem. 
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2 

However,  Q57  will  always  be  a  positive  quantity  unless  it 
is  defined  as 

Q§ 7  “  $57*  IQ57I  .  <«) 

in  which  case  the  sign  of  q|7  will  reflect  the  direction 

of  flow.  As  noted  previously,  K7  determines  the  sign  of 
the  final  term. 

2 

is  used  in  the  last  term  of  equation 
(40) ,  becaue  the  area  employed  is  always  that  of  the  smaller 
of  the  two  at  a  junction. 

Previously  undefined  variables  are 
Q57  -  <Q5  +  Q7)/2  (42) 

Le^7  =■  effective  length  from  (5)  to  (7)  . 
Effective  line  length  will  be  discussed  in  a  later  section. 

(3)  Node  Nine .  The  pressure  at  point  nine  has 
the  same  form  as  that  at  seven 

P9  -  P7-32u79Le79Q79/A7D^  -  1 . 14p79Q79/A7-p9K9Q9/2A7  (43) 

with 

Q79  *  Q79  I Q79 I  (44) 

for  the  reasons  discussed  in  the  previous  section. 

Flow  through  this  junction  is 

Q9  “  Q;q  "*■  Qei  (45) 

Qn  will  be  defined  with  respect  to  ram  flow  in  a  following 
segment  of  the  paper. 

Qei  "  CelP9  .  (46) 

The  coefficient  of  external  leakage  (Ce|)  is  the  allowable 
leak  rate  past  the  cylinder  tail  rods  to  the  atmosphere. 
Appendix  F  derives  the  values  used,  based  on  information 
contained  in  Ref.  25. 
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There  were  no  volume  terms  included  in 
equation  (45)  because  the  transition  from  the  end  of  the 
line  to  the  cylinder  entrance  was  considered  an  infinitesimal 
distance . 

b.  Secondary  Component  Effects 

It  is  readily  apparent  from  a  comparison  of 
Figure  (5)  with  Figure  (3)  that  the  simplified  system  omits 
numerous  valves,  fittings,  and  other  miscellaneous  elements 
which  are  installed  between  the  servovalve  and  ram.  Although 
these  "secondary  "  components  have  a  minor  influence  individually , 
when  combined  they  contribute  significantly  to  both  static 
and  dynamic  characteristics  of  the  system. 

The  approach  taken  in  this  development  was  to 
consider  them  in  three  categories:  branch  lines,  minor 
components,  and  major  elements. 

(1)  Branch  Lines .  Branch  lines  are  the  segments 
of  pipe  connected  to  the  main  line.  They  are  subjected  to 
system  pressure  but  do  not  ordinarily  allow  fluid  to  flow. 

Gage  lines,  hand  positioning  branches,  and  capped  sections 
were  included  in  Figure  (5) . 

These  lines  contribute  only  to  the  total 

system  volume. 

(2)  Minor  Components .  The  minor  components 
were  considered  to  be  any  elements  installed  in  the  system 
which  contributed  to  the  effective  length  and  volume  of 
the  system  but  were  not  specifically  listed  in  the  list  of 
major  components  [Ref.  26]. 
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The  elements  in  this  group  were  ball  valves, 
tees,  and  elbows.  Reference  27  was  used  as  the  source  for 
information  on  them  and  Appendix  G  explains  their  development 
for  use  in  the  overall  analysis  of  the  system. 

(3)  Major  Components.  The  deceleration  valves 
(decel)  and  the  blocking  valve  were  included  as  major  components. 
These  were  listed  in  Ref.  26  with  specific  pressure  and  flow 
characteristics.  The  volumes  employed  were  a  "best  estimate" 
based  on  diameter  and  the  author's  personal  experience.  The 
effective  lengths  were  as  determined  from  the  procedure  in 
Appendix  H. 

The  function  of  the  blocking  valve  is  to 
transfer  to  an  alternate  source  of  hydraulic  supply  (the 
Vital  System),  in  the  event  of  a  loss  of  main  header  pressure. 

As  pressures  in  the  two  systems  are  ordinarily  the  same 
(within  lOOpsi.  of  each  other)  and  the  valve  has  internal 
seals,  cross-port  leakage  is  negligible. 

The  deceleration  valve  is  a  cam  operated 
full-flow  valve  which  is  installed  to  prevent  over-travel 
of  the  ram.  As  the  ram  approaches  its  maximum  limit  of 
travel,  a  collar  on  the  tail  rod  mechanically  trips  the 
decel  can,  shutting  the  deceleration  valve.  This  permits  a 
fluid  cushion  to  arrest  the  ram  travel  rather  than  allowing 
it  to  a  "hard-stop"  (metal  -  to  -  metal)  at  maximum  velocity. 
Motion  in  the  reverse  direction  is  permitted  by  flow  through 
the  internal  check  valve  until  the  linkage  allows  the  main 
deceleration  valve  to  reopen. 
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(4)  Combined  Effects .  The  various  elements 
were  combined  in  tabular  form  for  inclusion  in  the  system 
analysis  (Tables  1  and  2) .  These  were  input  directly  into 
the  computer  program  to  solve  for  the  various  parameters 
as  discussed  above. 

3.  Cylinder  and  Piston 

The  final  element  in  the  system  is  the  cylinder  and 
piston  combination.  Here  the  power  is  transmitted  from  the 
fluid  to  the  stem  planes  by  the  piston  and  its  associated 
linkages . 

a.  Fluid  Dynamics 

The  pressure  and  flow  equations  assumed  a  form 
similar  to  those  of  the  lines. 


(1)  Flow.  The  flow  in  the  cylinder  was  described 


by 

<*11 

QU  +  ^11  +  Vll^ll/Bll 

(47) 

where 

*  Ci£(Pll"P12) 

(48) 

with 

P11 

-  / Pudt 

(49) 

Qu 

is  the  internal  leakage  from  the 

high 

pressure 

side  to  the 

low  pressure  side  of  the  piston. 

The 

determination  of  C^  is  shown  in  Appendix  F.  Flow  was 
defined  as  positive  from  node  (11)  to  node  (12);  thus, when 
**12  8reater  than  P^,  will  assume  a  negative  value 
as  required  by  the  line  equations . 
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TABLE  I  LINE  SEGMENT  DATA 


(50) 


'11 


AllXr 


^  is  the  effective  piston  area 


A11  -  *(Dll  -  Dr>/4 


(51) 


where  is  the  land  diameter  of  the  cylinder  and  Dr  is 
the  diameter  of  the  piston  tail-rod. 


X. 


r  is  the  velocity  of  the  ram. 
V11  =  A11L11 


(52) 


where 


L11  "  Lo  +  Xr 


(53) 


and 


L  =  %(L  -t  ) 
o  cy  p' 

Lcy  is  the  length  of  the  cylinder  t 


(54) 


is  the  thickness  of  the  piston.  Ram  position  was  defined 
as  positive  to  the  left  in  Figure  (5) ,  so  that  motion  to  the 
left  increases  while  movement  Of  the  piston  to  the  right 
decreases  V^.  For  this  reason,  it  was  necessary  to  define 


V10  =  L  -  X 
12  o  r 


(55) 


and  ultimately 


^12  =  ~^±l  ~  ^12  "  V12P12/612- 
Thus,  the  various  flows  in  the  system  would  automatically 

assume  the  proper  sign  to  reflect  the  direction  of  flow. 

(2)  Pressure .  The  pressure  equation  for  the 

cylinder  and  piston  arrangement  became 

P11  =  P9  '  32ullLllQll/AllDll  "  1*14  P11Q11/A11-  (57) 

There  was  no  entry  loss  term  because  it 

was  taken  into  account  with  Pg.  The  other  difference,  with 

regard  to  the  previous  line  loss  equations,  is  that  is  the 
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land  diameter  of  the  cylinder  while  is  the  effective 
area  of  the  piston.  The  effective  area  was  used  because 
Q/A  is  a  velocity  term  employed  in  defining  viscous  loss. 

On  the  other  hand,  D^  was  used  to  account  for  the  surface 
area  upon  which  the  fluid  is  "'dragging";  i.e.,  the  walls  of 
the  cylinder.  Since  the  tail  rod  is  traveling  at  essentially 
the  same  velocity  as  the  fluid,  it  contributes  no  losses 
due  to  viscosity. 

As  in  previous  sections,  it  was  necessary 

to  define 

Qll  “  Qn*lQiil  (58) 

to  account  for  flow  direction. 


b.  Piston  Dynamics 

Ram  motion  was  defined  in  terms  of  position  (Xr) , 
velocity  (Xr) ,  and  acceleration  (Xr) . 

Xr  *  /Xr dt  (59) 

Xr  =  JXrdt  (60) 

*  (AiiPn"Ai2P12‘moil^r"FL)/MT  (61) 

The  load  force  (FL) ,  as  described  at  the  outset, 

was  input  as  a  ramp  function  dependent  on  ram  position  and 

velocity. 

The  flow  of  oil  (raQi^)  term  was  included  to 
account  for  the  force  required  to  maintain  the  momentum  of 


the  oil.  It  was  determined  to  be 

®oil  =  p12^12 


for  &r>0 


(62) 


The  total  mass  of  the  system  (Hj.)  was  defined 

as  the  mass  of  the  ram  plus  the  mass  of  the  oil  from  the 

low  pressure  side  of  the  piston  to  the  return  tank,  such  that 

tank 

Hr  -  nL.  +  l  PiV.  (64) 

ram 

where  the  subscript  (^)  indicates  each  main  line  segment 
(flow  path)  from  the  ram  face  to  the  return  tank. 

C .  PROGRAM 

The  foregoing  equations  were  re-arranged  to  yield  the 
form  shown  in  the  attached  computer  program.  They  were 
modeled  in  a  Continuous  System  Modeling  Program  (CSMP  III) 
in  the  belief  that  it  was  the  best  modeling  method  available. 
However,  numerous  attempts  failed  to  accomplish  the  desired 
result. 

The  algebraic  equations  form  a  loop  for  which  no  solution 
was  found. 
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V.  SUMMARY 


A.  CONCLUSIONS 

It  is  felt  that  the  development  of  the  line  equations 
demonstrates  the  necessity  of  including  these  factors  in 
any  in-depth  analysis  of  the  submarine  hydraulic  system. 

It  is  also  the  author's  belief  that  the  approach  taken  in  this 
presentation  will  yield  worthwhile  results  for  any  system 
employing  extensive  quantities  of  piping. 

As  stated  in  the  introduction,  a  generalized  computer 
program  of  this  type  would  have  many  benefical  aspects  for 
both  the  submarine  designer  as  well  as  the  operator. 

B .  RECOMMENDATIONS 

The  first  recommendation  is  a  possible  means  of  clearing 
the  algebraic  loop  problem  currently  preventing  the  program 
from  executing.  This  would  require  clearing  all  of  the 
pressure  and  flow  equations  of  reference  to  any  variables 
except  for  themselves  or  ones  which  directly  invoke  an 
integral.  That  is 

Py  =  f  (  P9  79  » Pj  y  »  t  Py  >  ®tC  .  )  # 

However,  it  could  not  take  the  form  where  ?y  was  set  equal 
to  Pg  or  any  of  the  flow  terms  (Q^),  as  they  do  not  invoke 
an  integral  function. 

The  other  possibility  would  be  to  locate  another  program 
which  could  efficiently  handle  both  integral  as  well  as 
simultaneous  algebraic  equations. 
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APPENDIX  A 


T 


Determination  of  Servovalve  Volume 


No  values  were  obtainable  for  the  volume  of  the 
servovalve.  Therefore,  the  following  method  was  employed 
to  obtain  a  realistic  approximation. 

Using  nominal  dimensions  given  in  the  technical  manual 
[Ref.  20  land  the  knowledge  that  the  valve  is  made  of 
aluminum  alloy,  a  volume  to  weight  ratio  was  determined. 
Dimensions  18.5  x  8.0625  x  9.5  in. 

Specific  Gravity  of  alloy  approximately  2.8 
Given  weight  69  lb. 

Density  of  water  at  60°F.  62. 41b/ ft3 


Vbody 

Yalloy 

W  -  VY 


1417  in3 
0.101  lb/in3 
143.3  lb  for 


a  solid  block 


then 


wgiven/W 


0.48 


Thus  the  valve  body  weighs  approximately  half  of  what  a 
solid  block  would  weigh.  It  was  therefore  assumed  that 
the  volume  of  metal  in  the  valve  was  708in  .  Further, 
estimating  that  the  passages  for  the  pilot  system  occupy 
half  of  the  void  area  yielded  a  main  valve  volume  of  about 

3 

350  in  ,  and  since  this  volume  is  that  occupied  by  both  the 
supply  and  the  return  segments  of  the  valve,  half  was  assigned 
to  each. 

V,  -  Vo  -  175  in3 
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APPENDIX  B 


Null  Flow  Pressure  Coefficient  Determination 
The  allowable  null  flow  leakage  was  given  as  600  cc/min. 
at  1200  psi.  [Ref.  20  ].  This  was  converted  to  standard  form 

3 

of  (in  /sec) /psi.  as  follows 

Kco  *  600(cc.min)/1200  psi  x  1  in^/16.39cc  x  1  min/60  sec  • 
This  yielded 

K  =  5.0844  xlO-^  (in'Vsec) /psi. 

U 

and 

%Kcq  *  2.5422  xlO'^  (in^/sec) /psi  . 

APPENDIX  C 

Empirical  Valve  Leakage  Coefficient 
The  leakage  flow  across  a  valve  port  decreases  as 
the  valve  moves  to  block  the  port.  Since  a  plot  of  leakage 
flow  (Q^)  versus  valve  position  (Xv)  is  symmetrical,  with  an 
exponential  decay  as  shown  in  Fig.  C-l,  and  exact  data  was 
not  given  for  leakage  other  than  the  null  position  an 
empirical  expression  was  developed. 


Figure  C-l  Leakage  Flow  vs  Valve  Position 
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The  most  appropriate  expression  to  fit  this  curve  took 
the  form  of 

VW  e-^VXvmax)2 

where 

■  leakage  flow  at  any  valve  position 
Q^max*  maximum  leak  rate  *  %Kcq 
Xy max*  max^mum  valve  travel  from  null. 
Representative  values  were  then  chosen  such  that 

v<w  -  -02  VW  -  -25 

That  is,  at  25%  of  maximum  valve,  travel  leakage  through 
the  "blocked"  port  of  the  valve  decreases  to  2%  of  the  null 
position  leakage. 

Solving  for  A 

x  -  (‘"’WVI'IVW2 

A  -  62.6. 

In  reality  A  could  vary  from  near  zero  to  about  100. 
Realistically  it  would  have  a  value  between  50  and  100. 

A  -  65  was  chosen  as  the  value  to  employ  in  this 
presentation. 

The  form  used  was 

Q i  -  %Kcoe'X^Xv/Xvmax) 

with  the  determination  of  maximum  valve  travel  (X  _)  shown 

vm 

in  Appendix  D. 


APPENDIX  D 


Determination  of  Maximum  Valve  Travel 


Reference  shows  that 

Qmax  “  80  GPM  =■  80  x  231/60  in3/sec 
at 

T  *  115°F 

APtotal  "  1200  Psi 

APport  *  600  psi 
as  depicted  in  Figure  D-l. 


600  PSI 

Figure  D-l  Valve  Test  Flow  Path 

Using  the  orifice  equation 
Q  -  CdA(2Ap/p)^ 

with 

Cd  -  0.61 
p  -  8.442  xlO"5 
A  where  d^  *  1  in. 

*v  *  Xvm 

Q  “  Qm 
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such  that 


Xvm  - 

resulting  in  X„  ■  0.0426  in. 
vm 

This  value  was  used  in  spite  of  the  fact  that  it  seemed 
large  for  a  ''standard”  pilot  controlled  servovalve. 


r* 


APPENDIX  E 

Loss  Coefficient  Due  To  Abrupt  Area  Changes 

The  standard  form  of  the  pressure  loss  due  to  abrupt 
piping  diameter  change  is  given  as 
Hx  -  KV^g 

where 

is  head  loss,  expressed  in  feet  of  fluid 

V  is  fluid  velocity 

g  is  acceleration  due  to  gravity  • 

Head  loss  can  also  be  expressed  as 

El  -  Ap/a 

Y  *  Pg 

and 

y  is  weight  density 
while  velocity  is 

V  -  Q/A  -  4Q/ttD2 
combining  the  above  equations 
AP  -  8pKQ2/7r2D4  . 

K  is  defined  from  Figure  D-l 
K  -  %(1-D1/D2)  for  -K3 
K  -  -(1-(D1/D2)2)2  for  -Q. 

NOTE  is  always  the  smaller  of  the  two  diameters  being 
considered.  - -  ■  >~  Q+ 


Figure  E-l  Flow  at  an  Abrupt  Area  Change 
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APPENDIX  F 


Seal  Leakage  Coefficients 


Seal  leakage  was  based  on  permissable  values  listed  in 
Naval  Ship's  Technical  Manual, [Ref.  25.]  The  values  used  in 
this  presentation  were  for  °used'''  seals  with  maximum  allowable 
leakage  before  replacement  is  required.  This  was  felt  to 
be  reasonable  since  the  operating  temperature  of  the  oil 
is  normally  greater  than  the  test  temperature,  which  will 
produce  slightly  higher  lead  rates. 

The  leak  rates  were  given  as  one  milliliter  per  ten 
cycles  per  inch  of  rod  diameter  for  external  leakage,  and 
five  milliliter  per  five  minutes  per  inch  of  rod  diameter 
for  internal  leakage.  These  were  based  on  normal  operating 
pressures  and  a  fluid  temperature  of  100°F. 

Operating  pressures  were  assumed  to  be 
Pel  -  2000  psi 
Pu  *  1000  psi 

where  P£^  is  the  pressure  for  external  leakage  and  P^  is 
the  pressure  for  internal  leakage.  A  cycle  was  defined 
as  the  time  to  go  from  neutral  (zero  plane  angle)  to  full 
rise  to  full  dive  and  back  to  neutral  position.  For  the 
system  modeled,  stem  planes  travel  is  plus  or  minus  27°. 

Thus, a  full  cycle  would  be  four  times  this  value,  or  108°. 

The  normal  rate  of  travel  is  about  five  degrees  per  second. 
Thus,  108/5  -  21.6  sec/cycle  or  approximately 
3  cycles/min. 
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Since  the  tail  rod  used,  had  a  diameter  of  four  inches, 
the  allowable  external  leakage  was  four  milliliters  per  ten 


cycles . 
Then 


4ml  x  lcc  x  3_C2  .  i^cc/mt. 

10cy  1ml  1  min 


and 


1.2  cc 
min 


x 


3 

1  min  v  1  in  v 
60“i?c  TO?c 


1 

psi 


Ce£  *  6-1  (in3/sec) /psi 

Similarly  the  internal  leakage  for  a  land  diameter  of 
8.25  inches,  rounded  to  nine  inches,  was  nine  cubic 
centimeters  per  minute. 


And 

„  _  9 

CU  60  x  15739  x~TOO(J 

Ci£  *  9,15  xl0~6  (in3/sec) /psi 
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APPENDIX  G 


Minor  Component  Analysis 

The  values  used  in  the  derivation  of  the  minor  component 
effects,  considered  herein,  were  obtained  from  Crane  Technical 
Paper  No.  410,  [  Ref.  27] . 

The  volume  was  determined  to  be 
V  -  AL 

where  A  was  the  cross-sectional  area  of  the  pipe  entering 
or  leaving  the  element  and  L  was  a  value  determined  for 
each  type  of  element.  The  values  for  L  were  chosen  as 
Ball  valves  L  *  3  x  D 

Tees  L  *  2  x  D 

Elbows  L  *  1 . 6  x  D 

In  cases  where  the  element  has  differing  entry  and  exit 

diameters,  D  is  the  larger  of  the  two. 

The  values  for  the  effective  lengths  were 

Ball  valves  L  =  3  x  D 

e 

Flow  through  tees  Lg  *  20  x  D 

Branch  flow  tees  L  -  60  x  D 

e 

45°  elbows  L  ■  16  x  D 

e 

90°  elbows  L  ■  30  x  D 

e 

where  the  diameters  were  as  above.  An  additional  factor  of 
K  (appendix  E)  would  have  to  be  included  for  diameter 
differences  in  the  main  flow  path. 

The  value  of  effective  length  cited  above,  were  determined 
for  Reynolds  numbers  (Rg)  greater  than  1000.  Therefore, 


it  was  necessary  to  modify  them 


Since  Rg  =  4Q/(ttDv) 

and  L  L  R  /1000  for  R  <1000 

ed  ec  e  e 

with  ^ed*  desired  effective  length 

Lec“  ^e  calculated  effective  length  from  above 

then  Lfi  *  Lgc  (4Q/(-rrDv))xlO~3 

If  R  >  1000  then  L  =  L. 

e  -  e  ec 

so  that  for  4Q/ttDv>1000 

or  Qp/Dy >  785.4  L£  -  Lec 

and  for  Qp/Dy<785.4  Lg  -  1.273  xlO_3LecQp/Dy 
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APPENDIX  H 


Effective  Length  of  Major  Valves 


The  standard  head  loss  equation 
2 

H-^  =  Kv/2g  was  rearranged  in  the  form 
AP/pg  -  KQ2/2gA2 
This  yielded 

K  =  2aPA2/pQ2 

since 

A  =ttD2/4 
K  =rr2APD4/8pQ2 
and  since 

K  -  Le/D 
Le  ~tt2aPD5/8pQ2 

To  allow  Q  to  be  entered  in  the  form  presented  (GPM) , 
a  conversion  was  included. 

231  in3/gal  x  1  min/  60  sec 
so  that 

Q  (GPM)  x  231/60  =  Q  (in3/sec) 

and  Le  =  .0832  PD5/pQ2 
where 

D  =>  diameter  of  valve  (in) 

P  »  given  pressure  loss  (psi) 

Q  -  given  flow  (GPM) 

P  ■  density  of  fluid  at  given  temperature  (slug/in  ) 
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APPENDIX  I 


TO  EMPLOY: 

SUPPLY  LINES;  D1  IS  ALWAYS  THE  POINT  CLOSEST  TO  THE 
PUMP  -  THE  UPSTREAM  POINT 

RETURN  LINES;  D1  IS  ALWAYS  THE  POINT  FARTHEST  FROM 
THE  TANK  -  THE  UPSTREAM  POINT 

DUAL  DIRECTION  (REVERSING)  LINES;  D1  IS  THE  POINT 

NEAREST  THE  CONTROL  ELEMENT  - (VALVE) ,  VIEWED  AS  IF 
THE  LINE  ALWAYS  SUPPLIES  PRESSURE  TO  THE  POWER 
ELEMENT 

ENDMAC 

INITIAL 

CONSTANT  KCCCC=600 . , FT=1200 . , LAMB=65 . 

VALVE  LEAK  RATE  OF  MAIN  SPOOL 

NOTE:  LAMB  IS  STRICTLY  EMPIRICAL  AND  BASED  ON  AN 
AVERAGE  CURVE  OF  LEAK-RATE  VS  VALVE  DISPLACEMENT. 

CONSTANT  XVM=0 .042 5 , DLV=1 . 0 , CD=0 .61, Vl=l?5 . , V2=17  5 . 

constant  od5=i.9o,th5=o.25,l5=4S.o,ntb5=o. ,NTT5=0.,  ... 
nl45=o . , nl95=i • , nv5=o . 

CONSTANT  0D6=1 . 9 ,  TH6=0 . 25 ,  L6=l4-8 . 0 ,  NTB6=0 . ,  NTT6=0 . , . . . 

NL46=0 . , NL96=1 . , NV6=0 . 

CONSTANT  CD7=2 . 0 , TH7=0 . 27 , L7=228 . 0 , NTB7=2 . , NTT7=1 . , NL47=0 . , . . . 
NL97=6. ,NV7=1. 

CONSTANT  0D8=2 . 0 , TH8=0 . 27 , L8=2l6 . 0 , NTB8=3 . , NTT8=0 . , NL48=0 . , . . . 
NL98=5. ,NV8=1. 

CONSTANT  0D71=0 . 75, TH71=0 . 109 , L71=8 . 0 , NTB71=1 . , NTT71=C . .... 
NL471=0. ,NL971=0. ,NV71=0. 

CONSTANT  0D72=0.375,TH72=0.058,L72=120.0,NTB72=0. ,NTT72=0., . . . 
NL472=0 . , NL972=1 . , NV72=1 . 

CO NSTANT  0D7  3=1 . 5 , TH7  3=0 . 20  3 , L7  3=12 . 0 NTB7 3=0 . , NTT7 3=0 . , . . . 

NL47  3=0 . ,  NL97  3=0 . ,  NV7  3=0 . 

CONSTANT  0D8l=0 . 75 , TH8l=0 . 109 , L8l=8 . 0 , NTB8l=l . , NTT6l=0 . , . . . 
NL48l=0. ,NL98l=0. , nv8i=o. 

CONSTANT  0D82=0. 3750, TK82=0. 058, L82=120. 0, NTB82=0. ,NTT82=0. , .. . 
NL482. ,NL982=2. ,NV82=1. 

CONSTANT  0D8 3=1 . 5, TH8 3=0 . 20 3 , L8 3=12 . 0 , NTB8 3=0 . , NTT8 3=0 . 

NLl+8  3=0 . ,  NL98  3=0 . ,  NV8  3=0 . 

CONSTANT  PDR7=60 . 0 , QDR7=50 . 0 , RDR7=8 . 68E-05 , VDV7=60 . , DDV7=1 .094 
CONSTANT  PER7=50 .- ,  QBR7=50 . ,  RBR7=8 . 68E-05 ,  VBV7=20 . ,  DBV7=1 . 46 
CONSTANT  E=3.0E07,P0IS=0.3 
CONSTANT  B0=200000. 

CONSTANT  PR=50. ,PIN=3000. , FA=0 . 1 , 0DRE=2. 0 , THRE=0 . 148 , VR-2750 . 
CONSTANT  0ILCC=9 . 0 , CELCC=1 . 2, PCEL=2000 . , PCIL=1000 . 

CONSTANT  DLR=8 . 255 , DRR=4»0 , 0DR=11 . 5 , LCY=39 .875, TR=5 . 125 , MR=12 . 1 
CONSTANT  PI=3. 141593 
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KC  c=0 . 5*KC0  C  C/ ( PT*98  3 . 4 ) 

THE  VALUE  OF  (0.5)  IS  USED  BECAUSE  THE  TECH.  MANUAL 
LISTS  A  LEAK- RATE  FOR  THE  WHOLE  VALVE;  THE 
VALUE  SOUGHT  IS  THAT  FOR  HALF  THE  VALVE  (ONE  SET 
OF  PORTS). 

KCOCC  IS  TECH.  MANUAL  LEAK-RATE  IN  (CC/MIN). 

PT  IS  TEST  PRESSURE  USED  TO  ESTABLISH  THE  LEAK-RATE  (PSI). 

KCO  IS  LEAK-RATE  IN  ( ( IN**3/SEC)/PSI ) . 

AG=CD*PI*DLV*SQRT ( 2 . ) 

AG  IS  THE  PRODUCT  OF  THE  CONSTANT  COEFFICIENTS  USED  IN 
IN  THE  SERVO -VALVE  EQUATION. 

CFA=FA/1.4 

FA  IS  THE  FRACTION  OF  AIR  ENTRAPPED  IN  THE  SYSTEM, 

EXPRESSED  AS  A  PERCENTAGE  OF  THE  TOTAL  SYSTEM  VOLUME. 
CFA  IS  THE  COEFFICIENT  OF  THE  AIR  TERM  IN  BULK  MODULUS 
CALCULATIONS. 

k5p , k5m=snt  co ( di , d5 ) 

K6p , K6M=ENTC0 ( D2 , D6 ) 

K7P,K7M=ENTCO(D5,D7) 
k8p ,k8m=entco ( D6 , d8 ) 

K9P , K9M=ENT  CO ( D7 , Dll ) 

K10P,K10M=ENT  C0(D8,D12) 

THESE  COEFFICIENTS  ACCOUNT  FOR  EXIT  LOSSES  ONLY; 

ENTRY  LOSSES  ARE  CONSIDERED  AS  EXIT  LOSSES  FOR 
THE  NEXT  ELEMENT  UPSTREAM. 

vp5 , a5 , lep5 , d5la5 , bcs5=vale ( p i , e , po is , od5 , th5l5 , ntt5 , ntb5 , . . . 

nl45,nl95,nv5) 

VP6 , A6 , LEP6 , D6LA6 , BCS6=VALE ( PI , E , PO IS , 0D6 , TH& , L6 , NTT6 , NTB6 , . . . 

NL!|i>,NL96,NV6) 

VP7,A7,LEP7,D7,LA7,BCS7=VALE(PI,E,POIS,OD7,TH7,L7,NTT7,NTB7,. .. 

NL47»NL97,NV7) 

vp8,a8,lep8,d8,la8,b(s8=vale(pi,e,pois,od8,th8,l8,ntt8,ntb8, ... 

NLl+8 ,  NL98 ,  NV8l ) 

V&B1 , A7B1 , LE7B1 , D7B1 , LA7 Bl , BCS7 1=VALE ( PI , E , P0 IS , 0D7 1 , TH7 1 , . . . 

L71 , NTT71 , NTB71 , NL471 , KL971 , NV71 ) 
V7B2,A7B2,LE7B2,D7B2,LA7B2,BCS72=VALE(PI,E,P0IS,0D72,TH72, . . . 

L7  2 , NTT7  2 , NTB7  2 , NLU7  2 , NL9 7  2 , N V7  2 ) 
V7B3,A7B3»LE7B3,D7B3,LA7B3>BCS73=VALE9PI,E,POIS,OD73,TH73,.. . 

L7  3 ,  NTT7  3 ,  NTB7  3 ,  NLii-7  3 ,  NL97  3 ,  NV7  3 ) 
V8Bl,A8Bl,LE8Bl,D8Bl,LA8Bl,BCS8l=VALE(PI,E,POIS,OD8l.TH8l, . .. 

LSI ,  NTT81 ,  NTB81  ,NLi|8l ,  NL9ol ,  IIV8I ) 
V8B2,A8B2,LE8B2,D8B2,LA8B2,BCS82=VALE9PI,E,P0IS,0D82,TH82, . . . 

L82 , NTT8  2 , NTB82 , NL482 , NL982 , NVb  2 ) 
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v8b3,a8B3,LE8B3,d8b3,LAQb3,BCS83=VALE(PI.POIS,OD83,TH83,... 

l83,NTt83,NTB03,NL483,NL983,Nv83) 

B0I=1.0.B0 

BCII=2.0*(1.0+P0IS)/E 

BI1=B0I+BC1I 

BI2=BI1 

B15~BI1 

BI6=BI1 

BI7=BI1 

Bl8=BIl 

BI9=BI1 

BI10-BIX 

BIH=BI1  • 

BI12=BI1 

bcdv8=bii*vdv8 

BCDV7=BIl*VDV7 

bcbv8=bii«-vbvo 

BCBV7=BI1*VBV7 

bi57=boi+bcs5/v57 

Bl68=BOI+BCS6/V68 

BI79-BOI+ ( BCS7+BGS71+BCS72+BCS7  3+BCS74+BCDV7+BCBV7 ) /V79 
Bl8l=BOI+(BCS8+BCS8l+BCS82+BCS83+BCS8l4.+BCDV8+BCBV8)/V*l 

CQ57 , cq57s=linco ( a5 , d5 ,0d5) 

C058 , CQ66s=LINGO ( A6 , d6 , OD6 ) 

CQ79, CQ79S=LINGO ( A7 , D7 , OD7 ) 
cq8i, cq8is=linco ( a8 , d8 , od8 ) 

CQ11, CQllS=LINCO ( All , Dll , ODR) 

CQ12 , CQ12S=LINC0 ( A12 , D12 , ODR ) 

PlZ=(PIN-PR)/2,0 

P2Z=01Z 

P57Z=P1Z 

P68z=P1Z 

P79Z=P1Z 

P81Z=P1Z 

P11Z=P1Z 

P12Z=P1Z 

Q7Z=0.0 

q8z=o.o 

Q9Z=0  .0 

Q10Z=0.0 

Q11Z=0.0 

Q12Z=0.0 

Q57Z=0.0 

Q68Z=0.0 

Q79Z=0.0 

Q8lZ=0.0 

A11=PI*(D11«-D11-D11R*DHR)/U.0 

Al2=PI*(Dl2tfDl2-Dl2R*Dl2R)/4.0 
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D11=DLR 

D12=DLR 

DUR-DRR 

D12R*DRR 

a5s*a5*a5 

A6s=a6*a6 

A7S=A7*A7 

a8s=a8«-a8 

L57=La5 

A68=LA6 

L79-LA7 

L8l=LA8 

lep57*lep5 

LEP79=LEP7+LDV7+LBV7 

LEP68=LEP6 

LEP8 1=LEP8+LDV8+LBV 8 

V57=VP 5 
V68=VP6 

Y79=VP7+V7B1+V7B2+V7B3+VDV7+VBV7 

v8i=vp8+v8bi+v8b2+v8b3+vdv8+vbv8 

LDV7-SPECV ( PDR7 , QDR7 , DDV7 ) 
LDV8=LDV7 

LBV7=SPECV ( PBR7 , RER7 , DBV7 ) 

LBV8=LBV7 

VBV8=VBV7 

VDV8=VDV7 


DRE=ODRE- 2 . 0*THRE 

BIR=B0I+2 . 0-::- ( ( 1 .  +P0IS )*0DRE*0DRE+  ( 1 .  -P0IS ) *DRE*DRE) / . . . 
(E*(ODRE*ODRE-ORE*ORE) ) 

CEL=CELCC  /(PCEL*983.4) 

CIL=CILCC/  ( PCIL-”-98  3 .  U 

XRZ=0.0 

ZRDZ=0.0 

D1=3.*DLV 

D2=D1 

LC=0.5*(LCY-TR) 

DYNAMIC 
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«  PROCEDURE  SERVO- VALVE  *## 


PROCED  Ql,Q2,Q3,Q4*SERVO(AG,KCO,LAMB,PIN,PR,Pl,P2,... 

XV , XVM , RHOl , RH02 ,RHCR) 

CEX=KCO*EXP ( -LAMB* ( XV/XVM) **2 ) 

IP(XV)10,20,30 
10  Q1=CEX*(PIN-Pl) 

Q2=AG*SQRT { XV*XV* ( PIN- P2 ) /RH02 ) 

Q3=AG*SQRT ( ZV*ZV* ( ?1=PR)/RH0R) 

Q4-=CEX*(P2-PR) ' 

GO  TC  kO 

20  Ql=CEX* ( PIN-Pl) 

Q2=*CEX*(PIN-P2) 

Q3=CEX*(Pl-PR) 

Q4.=CEX*(P2-PR) 

GO  TO  4.0 

30  Ql=AG#XV*SQRT ( (PIN-PlJ/RHOl) 

Q2=CEX*( PIN-P2) 

Q3=CEX*(P1-PR) 

Q4=AG(XV(SQRT( (P2-PR)/RH0R) 

^0  CONTINUE 
ENDPRO 

**  PROCEDURE  KPRO  *** 

PROCEED  K5 , K6 , K7 ,K8 , K9 , K10=KPR0 ( Pi , P2 , P5 , F6 , P7 , P8 , P9 , P10 ,K5P , • . 
K5M,K6P,K6M,K7P,K7M,K8P,K8M,K9P,K9M,K10P,K10M) 
if(pi.ge.f5)go  to  100 
k5=k5m 
GO  TO  101 

100  k5=k5p 

101  IF(P2.GE.P6)G0  TO  102 
K6=K6M 

GO  TO  103 

102  K6=K6P 

103  IF  (P5.GE.P7) GO  TO  104. 

K7=K7M 

GO  TO  105 
104.  K7=K7P 

105  IF(P6.GE.P3)G0  TO  106 
K=K8M 

GO  TO  107 

106  K8-K8P 

107  IF(P7.GE.P9)G0  TO  108 
K9=K9M 

GO  TO  109 

108  K9=K9P 

109  IF(P8.GE.P10)G0  TO  110 
KlOsKlOM 

GO  TO  111 

110  K10-K10P 

111  CONTINUE 
ENDPRO 
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***  PROCEDURE  EFFECTIVE  LINE  LENGTH  *** 


PROCED  LE57=PR057 ( LEP57 , L57 , Q57 , D5 , RH057 , MU57 ) 
CLE57=q57*RH057/ (mu57*d5 ) 

IF( CLE57.LT. 785.i4-) GO  TO  200 

le57=l57+lep57 

GO  TO  201 

200  LE57=L57+1.273E-0>CLE57*LEP57 

201  CONTINUE 
ENDPRO 

PROCED  LE68apR068 ( LEP68 , L68 , q68 , D6RHO68 , MU68 ) 
CLE68=Q68RH068/(MU68*D6 ) 

IF(CLE68.LT.785.l4*)G0  TO  202 

LE68=l68+lEP68 

GO  TO  203 

202  LE68=L68+1 . 27  3E-0  3*CLE68*LEP68 

203  CONTINUE 
ENDPRO 

PROCED  LE79-PR079(LEP79,L79,Q79,D7,RH079,KU79) 

CLE79=Q7  9*RH07  9/  ( MU7  9*D7 ) 

IF (CLE79.LT. 785. 14.) GO  TO  204 
•  LE79=L79+LEP79 

GO  TO  205 

204  LE79=L79+1. 273E-03*CLE79*LEP79 

205  CONTINUE 
ENDPRO 

PROCED  LE81=PR081(LEP81,l81,q81,d8,RH061,MU81) 
CLE81=q81-::-rho81/(MU81*d8  ) 

IF(CLE8l.LT785.i4*)G0  TO  206 

LE8l=L8l+LEP8l 

GO  TO  207 

206  LE81=l81+1.273E-03*CLE81*LEP81 

2 07  CONTINUE 
ENDPRO 

***.  PROCEDURE  MASS  Op  OIL 

PROCED  MD0IL,MT=K0PR0(MR,XRD,RH011,RK012,V11,V12,V11D, 
V12D, RH057 , RH068 , RH079 , RH081 , RH01 , RH02 , VI , V2 , V57 , Vb 
V79, V8l, . . . 

RHOR,VR) 

IF(XRD.GE.O.)GO  TO  I4.OO 

N0IL=RH011'-*V11+RH079*V79+RH057*V57+RH0*V1+RH0Rv-VR 
GO  TO  401 

I4.00  MOIL=RHOl2*Vl2+RHO8l*v8l+RHO60*RHO2*V2+RHOR*VR 

MD0IL=RH012*V12D 
401  MT=MR+K0IL 

ENDPRO 
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P1D=(Q1-Q3-Q5)/V1*B1I) 

P2D= ( Q2-QU-Q5 ) / ( V2*B2I ) 

P57D= ( oS-Ql ) / ( B57l*v57 ) 

P68d= ( Q6-Q8 ) / ( b68i*v68 ) 

P79D=  ( Q7-Q9 )  /  ( B79I-3-V79 ) 
p81D=(q8-Q10)/(b81I-::-V61) 

PllD= (Qll-QlL-VllD) / ( V11I*V11 ) 

P12D=(Q12+QIL+V12D ) / ( Bl2I*V12) 

P5=2.0*p57-P7 

P7=2.0*P79-P9 

P9=(Q9-Q11)/CEL 

p6=2.0*p68-p8 

P8=2.0-*P81-P10 

PlO= ( Q10-Q12 ) / CEL 

***  LINE  AND  RAM  FLOWS  *** 

q5=a5#sqrt ( 2 . o* ( pi - p5 ) / ( rho  5*k5 ) ) 

q6=a6*-SQRT  ( 2 . 0*  ( P2-  P6 )  /  ( RHO  6*K6 ) ) 

THE  FOLLOWING  EQUATIONS  IN  IMPLICIT  FORM  WERE  INSERTED  IN 
AN  EFFORT  TO  CLEAR  THE  ALGEBRAIC  LOOP  PROBLEM.  IT  DID  NOT 
WORK. 

THE  NEXT  LOGICAL  STEP  WOULD  BE  TO  COMPACT  ALL  EQUATIONS 
SUCH  THAT  EACH  VARIABLE  IS  EQUATED  ONLY  TO  ITSELF  OR  A 
VARIABLE  WHICH  IS  DIRECTLY  COUPLED  TO  AN  INTEGRAL. 

Q7=IMPL(Q7Z,O.OOI,Q7CC) 

Q7C=*(  (P5-P7-CQ57S*RH057*ABS(o57(-0.5»RH07*K7Q7Q7/Q5S)  . . . 
/(CQ57*MU57-::-LE57)  )*2.-q5 

Q7CC=Q7C 

q8=impl(q8z,o.ooi,q8cc) 

q8c=( (P6-p8-cq68s*RH068#q68#ABS(q68)-0.5*RH08*k8*q8/a6s) . . . 

/ ( CQ68*MU66*LE68 ) ) *2 . -q6 

q8cc=qQc 

Q9=IMPL(Q9Z,0.001,Q9CC) 

Q9C=(  ( P7-P9-CQ79S*RH079»Q79*ABS  ( Q79 )  -0. 5*RH09-::-K9*Q9*Q9/A7S. . 
/  ( CQ7  9*HU7  9-”LE7  9 )  )*2 .  -Q7 

Q9CC=Q9C 

Q10=IMPL(Q10Z,0.001,Q10CC) 

QI0C=  ( (  p8-Pl0-CQ8lS*RK08l*Q8l*ABS ( q81  ) -0 . 5*RH0l0-::-K10*Ql0*Ql0 
/ a8  S ) / ( CQ8 1*mu8 1« LE81) )*2.-q8 

Q10CC=Q10C 

Q11=IKPL(QHZ,0.001,Q11C) 

QllC=(P9-Pll-CQllS-::RH0ll::Qlli:-ABS(Qll))/(CQll*KUllvLll) 

Q12=IMPL(Q127,0.001,012C) 

Q12C=  ( P10- P12-CQ125*RK012*Q12«-ABS  ( Q12 ) )  /  ( CQ12*MU12<:-L12 ) 
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Q57=IMPL(Q57Z»0.0001,q57CC) 

q57c=(p5-P7-cq57s*rho57*q57*abs(<^7)-o.5*rho7*k7*q7-*Q7/a5s)  . . . 

/ ( cq57*mu57*le57 ) 
q57cc*q57c 

Q68=IMPL ( q68z , 0 . 00 1 , q68  CC ) 

Q68c*(P6-P8-CQ68s#RH068«ABS(q68)-0.5*RH08*k8*q8*q8/a6s) . . . 

/ ( CQ68*MU68*LE68 ) 

Q68CC=Q68C 

Q79=IKPL(Q79Z,0.001,Q79CC) 

Q79C=(P7-P9-CQJ9S*RH079*Q79*ABS(Q79)-0.5*-RH09*K9*Q9*Q9/A7S) . . . 

/(CQ79*KU79*LE79) 

Q79CC=Q79C 

Q81=IKPL(q81Z,0.001,q81CC) 

Q8ic=  ( P8-P10-CQ81S*RII081*ABS ( Q8l )  -0 . 5*RH010*K10*Q10*Q10 . . . 

/  a8s  )  /  ( CQ8 1*MU8  l-::-LE8 1 ) 

q8icc=q8ic 

QIL=CIL»(P11-P12) 

BlI=CFA/Pl+BIl 

B2I=CFA/P2+BI2 

BRI=CFA/PR+BIR 

b5i=cfa/p5+bi5 

B6i=CFA/P6+BI6 

B7I=CFA/P7+BI7 

b8i=cfa/p8+bi8 

B91=CFA/P9+BI9 

B10I=CFA/P10+3I10 

BllI=CFA/PU+BIll 

B12I=CFA/P12+BI12 

b57I=cfa/p57+bi57 

B68l=CFA/P68+Bl68 

B79I=CFA/P79+BI79 

B8ll=CFA/P8l+Bl8l 

•s-K-a-R- DENSITY  VALUES 

RHOR=8.256E-05+8.l4.23E-05-::-PR*BRI 

R0Hl=8.256E-05+8.i4,23E-05*Pl*BlI 

RH02=8 . 256E-05+8  Jl2  3E-05j.:-P2*E2I 

RH05=8 . 256E-05+8 . U2 3E-0  5*p5*b5i 

RfI06=8.256E-05+3.1+23E-05-::-p6-::-E6l 

RHO7=8.256e-05+8.1+23E-05*P7«-B7I 

RH08=8.256E-05+8.ii.23E-05-::-p8*B8l 

RH09=8 . 256  E-05+8 .  1+2  3E-05---p8v-B9l 

RH0l0=8.256E-05+8.l4.23E-05*P10*B10l 

RH011=8 . 256E-05+8 .  it  23E-05*P11*E1 1 1 

RH0l2=8.256E-05+8.l|23E-05vPl2-u-Bl2I 

RHO57=8.256E-O5+3.U23E-05v-P57*E57I 

RH068*8 . 2$6e-05+8  .  2  3E-05*P68*b68i 

RH079=8 . 256E-05+8 . i*23 E-05*P7  r*-B79I 

RHO81=8.256e-O5+6.14.23E“05^P81«-e81I 
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VISCOSITY  EQUATIONS  *** 


MU57=fl. 198e-06*EXP ( P57*B57I ) 
HU68=8 . 198e-06*EXP ( p68*B68l ) 
MU79=8 . 198E-06*EXP ( P79*B79I ) 
MU81=8.198e-06*EXP(p81-»b81I) 
MU11*8 . 198e-06*EXP  ( P1K-B11I ) 
MU12=8 . 198E-06*EXP ( P12*B12I ) 
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